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^ ■ Abstract 

We study the effect of electromagnetic radiation on the condensate of a Bose gas. In 
an earlier paper we considered the problem for two simple models showing the cooper- 
ative effect between Bose-Einstein condensation and superradiance. In this paper we 
formalise the model suggested by Ketterle et al in which the Bose condensate parti- 
cles have a two level structure. We present a soluble microscopic Dicke type model 
describing a thermodynamically stable system. We find the equilibrium states of the 
system and compute the thermodynamic functions giving explicit formulae expressing 
the cooperative effect between Bose-Einstein condensation and superradiance. 
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1 Introduction 



The present paper is motivated by the recent experiments exhibiting a special coherent in- 
teraction between matter and light, which has been nicknamed "four- wave mixing" jlj. In 
these experiments boson atoms with an internal structure, condensed in a trap, are irradi- 
ated with light produced by an external laser beam. The structure of the atoms is usually 
represented by considering them as having two levels |2J 0] • A system of two level-atoms 
interacting with light is very reminiscent of the Dicke model jl]. Moreover an important 
feature of this model, namely superradiance has been observed in these experiments, where 
it is found that there is an enhancement of both Bose-Einstein Condensation (BEC) and 
light radiation (superradiance) due to the interaction. 

Recently various models pQ - [3] for this BEC-superradiance coupling were constructed and 
discussed in order to describe both equilibrium and non-equilibrium superradiance by con- 
densed atoms. It is interesting to note that as early as 1978 Girardeau pj| had already 
anticipated this phenomenon in the context of superfTuid helium and had discussed the pos- 
sible impact of the equilibrium superradiance on the thermodynamic properties of the latter. 
In our recent letter jB] we have considered two simple systems by which we modelled the 
coherent behaviour of the BEC atoms irradiated by a laser beam, showing rigorously that a 
weakened form of the "four-wave mixing" interaction enhances the superradiance and BEC 
as proposed by Ketterle et al [I] - [3]. 

The aim of the present paper is to consider a model which takes explicitly into account the 
internal structure of the boson atoms. In fact we assume that our bosons have an internal 
two-level structure of the type described by SU(2)-spm symmetry. Therefore the one-particle 
wave functions are of the form i[)<S>s where ip G L 2 (R U ) describing the spacial localization and 
sgC 2 describing the internal (spin) state. Only the condensate particles, i.e., the particles 
in the ground state are supposed to interact with the external field, and therefore only the 
ground state boson particles are given a different ground state energy parameterized by a 
separation level parameter e. If e is put equal to zero, it is as if we have just two different 
types of boson particles. The interaction turns out to be a second quantized version of the 
well known Dicke maser model. In our model we suppose that the recoil of the particles 
is negligible. The model is in fact a realization of the physical mechanism explained in [2]. 
For our model we study the equilibrium states in the infinite volume limit (thermodynamic 
limit) and compute the corresponding thermodynamic functions. We examine the presence of 
cooperation between the BEC condensate and superradiance as a function of the separation 
level parameter e. The existence of this phenomenon confirms the results obtained in jU] for 
a simpler model. It can be seen explicitly from the expressions for the occupation densities 
for the bosons and photons. Our results predict that with conventional BEC one obtains the 
same phenomenon of BEC-superradiance cooperation as is observed for trap experiments. 

We note that experimentally one can observe the photon recoil effect which, on light atoms, 
can be non-negligible pp. However in the present paper we consider the case when the photon 
momentum is very small so that the recoil effect can still be neglected. In a later publication 
we shall study another model in which the influence of recoil is included. 



2 The model and its equilibrium states 



We consider a system of two types of bosons of mass m enclosed in a cubic box A in v 
dimensions (A C W) with volume V, centered at the origin. 

As usual let A* = {2nk/V^ u \k G Z u } be the dual space of A used to formulate the model 
with periodic boundary conditions. For k G A*, a — ±, a£ CT and the usual boson 

creation and annihilation operators of the two types of bosons satisfying the commutation 
relations: 

a *k',a'\ — fik,k'8<r,<r'- (2-1) 

The kinetic energy of the system is given by 

T a = Y1 Yl e ( k ) a l,* a k,a + £(a* + a 0)+ - ao,- a o,-) , (2.2) 
cr=± feeA*,fc^o 

where e > and e(k) = \\ k\\ 2 /2m. Note that the two k = mode bosons (the ground state for 
non-interacting bosons) have a supplementary internal energy, a spin-state energy, making 
the internal structure of the bosons explicit. On the other hand the excited bosons k ^ 
are not distinguished by their internal energy, but it is straightforward to make them also 
distinguished. The reader will able to see that that our arguments cover also the situation, 
when the single particle boson spectrum is presented by two branches: e a (k) := e(k) + as 
for two internal states os bosons. 

We represent the external one mode laser field by a single mode boson field with creation 
and annihilation operators b, b* satisfying [b,b*] = 1. As we indicated in the introduction 
here we consider the case, when the photon momentum is very small so that the recoil effect 
is negligible. In this approximation we can take k = and then 



1 



dxb(x), (2.3) 



A 



where b(x),x G W stands for the local (annihilation) photon field. As suggested in the 
introduction we define our model Hamiltonian as 

H A = T A + U A (2.4) 

where 

U A = -^=(a* 0+ a _b + a 0+ a* „b*) +Qb*b+ ^-N 2 A (2.5) 
2v V 



and 

N A =J2 N k, N k = (N kt+ + N k ,_) , N Ka = al a a k , a 
fceA* 

are respectively total boson number operator, the /c-boson number operator and the boson 
number operator for momentum k and type a. 

Furthermore SI > is the laser frequency and g is the coupling constant of the interaction 
between the bosons and the external field. Note that without loss of generality we can 
take g to be positive as we can always incorporate the argument of g into b by a gauge 
transformation. 



Notice in (|2.5jl the presence of the mean-field repulsive particle interaction with a positive 
coupling constant A > 0. This term is essential in order to obtain a model describing a 
thermodynamically stable system, i.e. ensuring the right thermodynamic behaviour. Indeed 
one can check by considering the interaction U A in (J2.5j) . that 



On the basis of the trivial inequality 4ab < (a + b) 2 , the lower bound of (|2.fi|) is bounded 
from below, if A > g 2 /8Q, that is if the stabilizing repulsive interaction coupling constant A 
is large with respect to the coupling constant g or if the laser frequency Q is large enough. 
Therefore we assume that A > g 2 /8Q is satisfied for the model (|2.4|) . The reader will see all 
along in the explicit analysis of the model below, the importance of this stabilizing condition. 
We note that, so far, neither the coherent recoil model, nor the "four- wave mixing" model 
nor Girardeau's model are thermodynamically stable, although Girardeau in [S] has stressed 
the importance of this stabilization. The models in jB] are stable because of the linearity of 
the interaction. 

In the present paper we study the equilibrium states of the model (|2.4jl in the grand-canonical 
ensemble and therefore we shall work with the Hamiltonian 

H A (fi) = H A -nN A (2.7) 

where \x is the chemical potential. Specifically our objective is to identify the infinite volume 
equilibrium states corresponding to the Hamiltonian (|2.7|) for a system of three different 
types of bosons. One way of achieving this goal is through the basic variational principle of 
statistical mechanics. Before starting to do this we prefer to reformulate the model with the 
purpose of showing that our model (|2.4|) is nothing but a second quantized bosonic form of 
the Dicke model and hence it realizes the ideas proposed in [2j and [B]. 

We have a system of atoms with internal states a — ±. Dicke regarded the two- level atom 
as a spin-1/2 system. This is what we shall also do and therefore we start from a two- 
dimensional representation of the Pauli matrices generating the Lie algebra of 577(2), given 
by 

'o i\ _ A) o\ 3 /io 



a -(0 oj' ff "=(i 0)' ° -{0 -1) (2 - 8) 

and the basis vectors {e+ = (1, 0), e_ = (0, 1)} of C 2 diagonalizing a 3 . 

The one-particle space of bosons is Ti, = L 2 (W) <g> C 2 . Let / <g> s be an element of 7i, then 

a*(f ® s) is the creation operator of a boson particle with state vector / ® s. In particular 
one can make the following identifications. 

4 j± = a*(/ fc <g)e ± ) (2.9) 

where for k G A*, f k is the plane wave function 

f k (x) = ^=e lkx , xeW. (2.10) 
yV 



In particular we have 

a* j± = a*(/o®e ± ). (2.11) 

For any <f> G TC, the creation and annihilation operators a*(4>) and a(4>) are linearly defined 
on arbitrary n-particle subspaces of Fock space TiTi): 

a*(0)sym(0! ® 2 ® . . . <g> <p n ) = (n + l)^sym(0 ® <j) X ® . . . ® <j) n ) (2.12) 

and 

n 

a(0)sym(0i <g> 2 <8> • • • <E> <fi n ) = ^(0, </v)h sym(0i <g> . . . <g> r <g> . . . <g> <j) n ), (2.13) 

r=l 

where sym denotes symmetrization, (•, is the scalar product in 7i, and r means that r 
is omitted. 

Applying these definitions for ± on the n-particle A; = mode states and using the identity 
<j + s = (e_,s)c2 e + we obtain 

ao+a - sym((/ ® si) ® (Jo ® s 2 ) ® • • • <8> (/o ® s n )) 

n 

= ^ + sym((/ ® si) <8> (/o ® s 2 ) ® ■ ■ ■ ® (Jo ® s n )) (2.14) 

r=l 

where 

vt(fo®si)(g)(fo(S)s 2 )(g). . .® (/ ®s„) = (/ ®si)®(/o®S2)®. • .®(/o®^ + s r )®. . .®(/o®s n ). 

(2.15) 

The fc = mode kinetic energy term can be treated similarly. Hence, on the n-particle k = 
mode states the sum of the k = kinetic-energy term ()2.2|) and the interaction term with 
the laser field ()2.5|) takes the form 

n n 

which coincides with the Dicke maser model. This proves that the model (|2.4j) (or (|2.7jl ) 
realizes the suggestions of 3], namely that it is nothing but a second quantized bosonized 
form of the Dicke maser model. 

So far we have discussed the structure of our model. The rest of the section is devoted to 
the technical preparation of the basic variational principle of statistical mechanics applied 
to our model (|2.7|) . 

The variational principle states that if S is the set of the extremal translation invariant 
states and / is the free energy density defined on S by 

f(u) = Jim oo(HM/V) - (l/0)S(u) (2.17) 

where S(u) is the entropy density of the state u>, then a state ujp G S satisfying 

/(^) = inf/H (2.18) 



is an equilibrium state of (|2.7jl at inverse temperature (3. 



The Hamiltonian ()2.7|) is not quadratic in the creation and annihilation operators, and 
therefore cannot be diagonalised by a standard symplectic or Bogoliubov transformation 
and thus, on this basis, one is tempted to conclude at first sight that the model is not 
soluble. However on closer inspection we find that we can write (|2.5jl in the form 

Ua g f / a* 0+ \ ( a - \ ( b \ ( a + \ / aj. 



v 2 lAvW \VvJ \VvJ \VvJ \WJ \vv 

b* \ ( b \ A /iV A x 2 



+" ' ~ 7 v) Wf) + 2 I ' ' < 219) 



so that all the terms are space averages. We have 
^= = i J^dxa±(x), ^| = ^|(iifl^) and ^ = i J^dxa* a (x) a{x) a . (2.20) 

and by virtue of (|2.3p . b* /y/V and bj\/V are clearly also space averages. Without going into 
all the mathematical details, the reason why space averages are such a simplifying feature is 
that they tend weakly to a multiples of the identity operator For example if u is a space 
homogeneous extremal (mixing) state then for all local observables, A and B one has 

lim u ( A — / dxa*(x) a(x)B ) = u(AB) lim uj ( — / dxa* (x) a a (x) ) 
V-*cc \ V J A J ; v-*oc \V J A 7 

= uj(AB)u(a*(0)a a (0)), (2.21) 
so that N A /V tends weakly to Y1 (J =± UJ ( a *(0)o- a o-(0))- Similarly 

lim ^± = w(o±(0)), and lim -?L = w(6(0)). (2.22) 



Thus if u; G S, then the contribution of the term (|2.5jl to the energy density in (j2.17j) yields 
^^V^ = 2 (^K(°)V( a -(0)V(&(0))+^(a + (0)V(a70)V(6*(0))} 

+fiK6(0))| 2 + ^g^«(0)a (J (0))j , (2.23) 

The result follows readily from (j2.21|) with A and B a multiple of the identity. We can 
therefore conclude that in the study of the equilibrium states of (|2.4J) or (|2.7|) . we can 
limit ourselves to searching for solutions u which are product states on the tensor product 
canonical commutation relations algebra (CCR) of the three different kinds of particles, 
namely on 

A := A+ <8> A- ® B, (2.24) 
where A± is the C* algebra generated by the Weyl operators: 

W±(/) := exp |z -j= 



for all / G L 2 (W) n L\R U ), and B by the Weyl operators: 

W b [f) : = exp jz — 1 . 

The above discussion makes it clear that we find the equilibrium states of our model amongst 
the states which are determined completely by their one-point and two-point functions, that 
is, among the set of extremal space invariant quasi-free states [7] on the respective CCR- 
algebras. This is a consequence of the fact that if u G S, the set of states on A, and Cj G S 
is a quasi-free state with the same one-point and two-point functions as u, then it follows 
by Klein's inequality [Sj that 

S(Q) > S(uj). (2.25) 

Therefore since our energy density involves only the one-point and two-point functions, if 
Sqf is the set of quasi- free state on A, then 



and consequently 



inf f(u) > inf f(u) (2.26) 



f{u p )= inf f(u). (2.27) 



We denote the set of quasi-free states on A a by u a determined by the constants a a and the 
non- negative operators A a on L 2 (IR I/ ) and satisfying 

= exp [W2 mt(a a (lj)) - i||/|| 2 - i</, A a f)^j (2.28) 

for all / G L 2 (W) n L 1 (R"), see 0. 

Note that the states io G are completely determined by the one-point function 

uj a (a a (f)) = aMA) (2-29) 

and the two-point function 

"MUM)) = (9, Kf) + \a a \ 2 (l,f)(g,l) (2.30) 
for all f,ge L 2 (W) H L^W). 

On B we consider the extremal invariant state, which is determined by one constant a^, see 




,(W b (f)) = exp (iV2 mt (a b (l, /» - I||/|| a ^ (2.31) 



Its one- and two-point functions are 

u b (b{f)) = a h (f,l) (2-32) 

and 

u b (b*(f)b(g)) = \a b \ 2 (lj)(g,l). (2.33) 

Note that this one-mode coherent state depends only on the k = mode. It is possible 
also to consider a more general quasi-free states of the form (12.31)) to take into account 



other photons modes. However since only the k = mode interacts with the bosons, this is 
unnecessary. 

Thus the candidates for the equilibrium states are among the set, Sp, of products of quasi- 
free states, i.e. states of the form 

uj = u + ® u_ (g> uj h . (2.34) 

They are completely parameterized by the set of parameters: a±, a\, G C and the integral 
operators A± on L 2 (M. U ): 



[ A ± {x-y)f{y)d»y, xeW (2.35) 



If A± is the Fourier transform of A±, then A±(k) > 0, expressing the positivity of the states 
u±. 

The variational principle (j2.18j) is now reduced to 

f(u p ) = inf f(u). (2.36) 

The entropy density for states in Sp is explicitly given by, see 0: 

S(u)=S(u+)+S(uJ), (2.37) 
where S(u)b) = because only one photon mode is taken into account. Here 

S(uj ± ) =J {(l + A ± (k)^j In (l + A ± (A:)) - A ± (k) \nA±(k)} 10-. (2.38) 

A straightforward computation yields: 



lim u(H\(n)/V) = —(/j, — e)|a + | 2 — (// + e)\a. 

V~*oo 

+ / {e(k) - (A+(k) + A-(kj) 



d v k 



(2nY 



+ ^ (a + a_ab + a + a_a{,) + Q \ab\ 2 



Note that the pressure P{n) of the system ()2.7|) . as a function of the chemical potential /x, 
is related to the grand-canonical free-energy density by 

P{fi) = -f{up) = -Mf{Lj). (2.40) 



3 Variational Solutions 



In this section we give a systematic derivation of the equilibrium states for our model as 
well as explicit expression for the corresponding grand-canonical pressure. To this end we 
solve the variational principal ()2.36|) . and we start by substituting ()2.37|) - (j2.39J) into ()2.36|) to 



obtain an expression for the functional f(oo) in terms of the variational parameters ai±, a& G 
C and A±(k): 

We find that there are two critical chemical potentials pP (e) and pf\e), Pc(s) < A*c (e)- 

For p < He (e), the two cr = ± Bose gases behave like two mean field Bose gases with no 
BEC and they do not interact with the external 6-boson laser field, in which there is no 
condensation either. 

For p > p^\e), there is BEC for the two a = ± Bose gases and for the external boson laser 
field (superradiance). 

When pi (e) < p¥\e), for p i c\e) < p < pf\e) there is BEC only for the a = — Bose gas. 

First we remark that we can take a±, ct^ real after a suitable gauge transformation on 
the boson creation and annihilation operators a ± and b, see f|2.29|) and (|2.32|) . Note that 
the squares of these parameters are in fact the condensate densities of the corresponding 
boson modes. For notational convenience we introduce the particle density for an arbitrary 
quasi-free state to of the form u + £g> oj- ® uit,, 

P := [ (Mk) + £± + K| 2 + M 2 = lim (3.1) 

J RV V / {2tt) u V^oo V 

that is, p is the density of cr = ± particles, excluding the 6-particles. We get the following 

Euler-Lagrange equations for the variational principle (|2.3(jj) : 

Take 

a + , cv_, a;, G M. (3.2) 

(i) Differentiation of f(uS) with respect to a + gives: 

2 (Xp + e — p) a + + ga-Ctb = 0, (3.3) 

(ii) differentiation with respect to a_: 

2 (Xp — e — p) cv_ + ga + ab = 0, (3.4) 

(iii) differentiation with respect to a^, 

2Qctb + ga + a_ = 0. (3.5) 

(iv) and finally differentiating with respect to A + and A_ yields: 

A + {k) = A_{k) = emk) _l +Xp) _ r (3.6) 

Note that the last equation implies that Xp — p > 0, since the A±(k) are positive. Moreover, 
the correlation inequality (see e.g. [TUj ) 

u([A*,[H A (p),A}})>0 (3.7) 

for all observables A, applied here with A = aj$_, implies that Xp — p > e > 0. Substituting 
(B into (jSU) we get 

p = \a + \ 2 + \a_\ 2 + 2p (p - Xp) (3.8) 



where 

pM := L e^W*)-M) - 1 (L)* (3 ' 9) 

is the density of the /ree Bose gas at chemical potential p. Recall that po(p < 0) < oo an d 
that p (p = 0) < 00 for > 2. 

Solving (J3.3j) . (J3.4j) and (J3.5j) we have to distinguish three cases: 
Case 1: a+ = a_ = cn& = 0. 

Substituting zero for a + and a_ into ()3.8j) we get the standard equation for density of the 
mean-field interacting bosons 

p = 2p (p-Ap), (3.10) 

see e.g. [O]. By virtue of the stability condition Ap — p > e > we see that this equation 
has no solution for p > pi(e) := 2Ap (— e) — e, while if p < pi(e) it has a unique solution 
p = pi(p). (See Figure [H where x = Ap — p so that 2 > e and ()3.10|) becomes p = 
2Apo(— x) — x). Putting this value of p into ([3.6)1 . we determine A±. Substituting these and 
a + = a_ = a.b = in the expressions (|2.31j) . (J2.33)) and (J2.33)) . for u± and u>b respectively, 
we find a solution, 0Up , of the Euler-Lagrange equations for the variational principle (12.36)1 
for p < pi(e). From ()3.8|) we are able to compute the free energy density for the state u4 • 



/(«£>,) = -2 Po (p - A Pl (p)) - ^Ap 2 (p) (3.1 I ! 



2 

where Po(a*) is the pressure of the free Bose gas: 

pM--=~ [ ^ (1 - e-W*)-")) (3.12) 

Case 2: a + , a_ and are non-zero. 
We obtain from (jH3J), fl32J and (J3~U) that 

2y/^(Ap- g -p) 2^(Ap + e-p) 2 V / (Ap-p) 2 -e 2 

a + = — - , a_ = — - , «t = . (3.13) 

9 9 9 

From these we see that in this case BEC is indeed present. Again substituting these values 
for a+, a_ into f|3.8|) we get 

p=^(Ap-p) + 2p (p-Ap). (3.14) 

g 2 

Note that the first term corresponds to the condensate density. Let 77 := (8QX/g 2 — 1). From 
the thermodynamic stability condition for (Section 2) we know that 77 > 0. Then equation 
(13.14)1 has a unique solution p = p2(p) for p > ^(e) := 2Apo(— e) + r/e. Substituting 
this value of p into 1)3.13)1 and (J3.6j) we obtain all the parameters a + , a_, and A± and 

(2) 

consequently we get another solution, Up , of the Euler-Lagrange equations. 

(2) 

The free energy density for the state ujp ' can again be computed: 

f{Jp 2 l) = -2 Po (p - Ap 2 (p)) - \\p\iti) + ^(Ap 2 (p) - p) 2 - (3.15) 



Denote by xq the unique solution of equation 2Xp' (— x) = r] corresponding to the minimum of 
the function 2Xp (— x) +r)x, and let p = 2Xp (—x ) + r]x . For p < p < ^{e) the equation 
f!3.14|) has two solutions p = p2{p) and p = f)2(p), P2(p) > P2i.i1)- The corresponding states 
u;^ and tDg can be found as above. The free energy density for the state ^\ is as in 
(j3.15|) and for Cj^ it is the same with P2(p) replaced by p 2 (/i). 



Case 3: 7^ and a + = otb = 0. 

From ()3.3|) . ()3.4|) and (|3.5|) one can see that this is possible only if 

P=—^~ ' ( 3 - 16 ) 

corresponding to the boundary x = e of the stability domain, see Figure ^ The equation 
()3.8|) then requires that p > pi(e) and gives 



a + = ^j^-2p (-e). (3.17) 

This case corresponds to yet another solution of the Euler-Lagrange equations, uj^ , whose 
free energy density is given by: 

f(J^) = -2p (-e) - (3.18) 

We see from above that for certain values of p there are several solutions of the Euler- 
Lagrange equations. Since these equations determine only the stationary points of the free 
energy functional, if there is more than one such point, in order to obtain the equilibrium 
state for a fixed p we have to decide which of the solutions, has the lowest grand-canonical 
free-energy density. 

To proceed with explicit analysis of solutions of the Euler-Lagrange equations it is easier 
to work with the variable x = Xp — p rather than p. Also in the grand-canonical ensemble 
it is more usual to use the pressure instead of the free energy density. These allow to find 
the grand-canonical pressure as a function of its natural variable, the chemical potential. In 
terms of x and 77 the equations (|3.10J) . ()3.14|) and (|3.16J) become: 

2Xp (—x)—x = p for p<pi(e), (3.19) 

2Ap (— x) + r/x = p for p > po (3.20) 

and 

x = e for p > pi{e). (3.21) 



We consider first the case e = 0. Then p±(0) = P2^0) = 2Xp c , where p c := Po(0). So, in 
this case the lower critical dimensionality the same as for the free (or mean-field) Bose-gas: 
v = 2. The equations (jSHHI), and (EQ) become : 

2Ap (-a;) - x = p for p< 2Xp C} (3.22) 



2Xp (— x) + r/x = p for p > po 



(3.23) 



and 

x = for p, > 2Xp c . (3.24) 

In Figured] we have drawn y = 2Xp (—x) — x and y = 2Xp (—x) + r\x. Recall that x is the 
unique solution of 2Xp' (— x) = r\ and po = 2Xp (—xo) + r/xo- It is easy to see that: 

1. For p < p , ()3.22|) has a unique solution Xi(p) while (I3.23J) does not have a solution. 

2. In the region p < p < 2Xp c , ()3.22|) has a unique solution Xi(p) while (j3.23|) has two 
solutions x 2 (p) and x 2 (/i), x 2 (p) > X2(p)- 

3. Finally for p > 2Xp c , ()3.22|) has no solution while ()3.23|) has a unique solutions x 2 (p) 
and we also have to consider the solution (j3.24j) . x = 0. 




y = 2Xp(-x)-x 



Let 



and 



Figure 1: Solution of the density equation 

(x + p) 2 



P 1 (x, t i) = 2p (-x) + 



2X 



P 2 (x,p) = 2p (-x) + 
Then the situation is as follows: 



{(x + p) 2 - (77 + l)x 2 } 
2X 



(3.25) 
(3.26) 



1. For /x < fi , the solution of the variational problem ()2.36|) is ui and the corresponding 
pressure P(/x) := -/(<«$],) = Pi (x x (p) , fi) ■ 

2. For /Xo < /x < 2Xp c , the solution of the variational problem is the state out of , 

(2) ~ (2) 

tjp and oj» which minimizes the free energy density or equivalently maximizes the 
pressure. The pressures for these states are Pi(xi(n), /x), P 2 (a; 2 (/x), /x) and P 2 (:r 2 (/x), /x) 
respectively. 

3. For /x > 2Ap c , the two candidates for the solution of the variational problem (j2.3fij) are 
u)^ and ■ The pressures for these states are P 2 (a; 2 (/x), /x) and P^{jj) : = —f(^f\) = 
2po(0)+/i 2 /2A. 




Figure 2: The pressure 



In Figure El we have sketched Pi(xi(p), /x), P 2 (x 2 (/x), /x), P 2 (x 2 (/x), /x) and P 3 (/x ). One can 
check that Pi(a;i(/x), /x) and P 2 (a; 2 (/x), /x) are convex in /x. One also has 



dP 2 (x 2 (//),//) X 2 (/i) + /! dP 2 (x 2 (fi),fi) X 2 (fi) + fi 



and 



d/x A d/x 

dPi(xi(/i), //) 



A 



d/x A 

Therefore since for /x < /x < 2Ap c , x 2 (/x) > x 2 (/x) > Xi(/x), in this interval we have 

dP 2 (x 2 (/x),/x) ^ dP 2 (x 2 (/x),/x) ^ dPi(xi(/x),/x) 



d/x 



d/x 



d/x 



(3.27) 



(3.28) 



(3.29) 



As P 2 (x 2 (p ), A*o) = P2(x 2 (p ), Po), it follows from ()3.29j) that P 2 (x 2 (p), p) > P 2 (x 2 (p),p) 
for p < p < 2Ap c . Now 

P 1 (x 1 (2Ap c ),2A Pc ) = P 2 (x 2 (2Ap c ),2Ap c ) = 2p (0) + 2Xp 2 c (3.30) 

and consequently P 2 (x 2 (2Ap c ), 2Ap c ) >Pi(xi(2Ap c ), 2Ap c ). Also if P 2 (x 2 (po), po) were greater 
than Pi(xi(po), po), then ()3.29|) would imply that P 2 (x 2 (2Ap c ), 2Ap c ) > Pi(xi(2Ap c ), 2Ap c ) 
contradicting ()3.30|) . Thus we must have 

PK^Po),/^) = P^^Po),/^) < P 1 (x 1 (p ),p ). (3.31) 

Therefore there exists a unique p c satisfying p < p c < 2Ap c such that P 2 (x 2 (p c ), p c ) = 
Pi(xi(p c ),p c ). 

Finally we consider p > 2Ap c . We have P 3 (2Ap c ) = Pi(x 2 (2Ap c ), 2Ap c ) < Pi(x 2 (2Ap c ), 2Ap c ) 
and 

dPM = P < (x 2 (p) +p) = dP 2 {x 2 {p),p) ^ 32 ^ 

dp A A <ip 

Therefore P 2 (x 2 (p),p) > P 3 (p) for p > 2Ap c . 

Summarizing: There exists a unique critical chemical potential p c such that 

1. For p < p c , the solution of the variational problem (j2.36|) is oj^ . For oj^ , a + = 
a_ = ctb = 0, i.e. the two ± Bose gases behave like two mean field Bose gases with 
no BEC and do not interact with the external 6-bosons which do not condense either. 
The corresponding pressure is P(p) = Pi(xi(p),p). 

( 2 \ 

2. For p > p c , the solution of the variational problem is u>p'. For this state 

2y/nx 2 (fi) 2x 2 (p) 
a + = a_ = — 1 and a& = , (3.33) 
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i.e. there is BEC for the two ± Bose gases and for the external bosons laser field (su- 
perradiance) . Moreover the condensation of the ± bosons is enhanced by the presence 
of the laser field (6-bosons), known as the equilibrium BEC superradiance The 
pressure for the system is P(p) = P 2 (x 2 (p),p). 

We now return to the case e > 0. We have to redefine P 2 and P 3 but Pi remains unchanged: 



and 

Note that 
and 



P 2 (x,p) = 2p (-x) + {( -^ )2 -^ +1) - 2} + ^ (3.34) 

P 3 (p)=2po(-^) + ^£. (3.35) 
P 1 (xi(pi(e)),p 1 (e)) = P 3 (pi(£)) (3.36) 



P 2 (x 2 (p 2 (s) ) , p 2 (e) ) = P 3 (p 2 (e) ) for e < x 

P 2 (x 2 (p 2 (f:)),p 2 (£)) = P 3 (p 2 (e)) for e > x . (3.37) 



Also by again considering the derivatives 



P 2 (x 2 (/i),fi)> P 3 (n) (3.38) 

for fi > ^t 2 (e) and as before P 2 (x 2 (p), fi) > P 2 (x 2 (p) , p) for fi < ft < 2\p (—e) in the region 
where it applies. 

(a) The simplest case to consider is when e > xq. In this case for ji < fii(s) only (j3.19|) has 
a solution £i(/i), for ii\(e) < \i < Ai 2 (e) only (j3.21|) is satisfied i.e. x = e and for ji > p 2 (e) 
only ()3.2()j) has a solution x 2 (fi). Thus the states are \ c^a and as fi increases. This 
means that as we increase \i the system goes from no BEC, to BEC for the a = — bosons 
only, to BEC for both species and superradiance. 

(b) When e < xq we have to consider two cases, fii(s) < po < ^{s) and /io < pi(s). 

In the first case we can use the same arguments as for e = to show that Ps(fio) > 
P 2 (x 2 (fio), (jl ) and P 3 (/i 2 (e)) < P2(x 2 (fj, 2 (e)) , fj, 2 (e)) . This implies that there exists fi c (s) 
between /i and fi 2 (e) such that P 3 (fi c (e)) = P 2 (x 2 (fj, c (e)) , fi c (e)) . Thus at n c (e) the state 

(3) (2) 

changes from to lu» . This means that the situation is the same as for e > x except 
that the changes of state occur at Hx{e) and at fJ, c (e). 

For £ < xo and /io < the same argument applies. However we did not determine 

on which side of the value of ix c {e) lies. Thus we know that there is no BEC for 

H < jiQ and there is BEC for both types of bosons for ji > fi c (e), but we do not know if the 
intermediate phase with BEC for — bosons only is present. 



4 Conclusion: Equilibrium BEC Superradiance 

The above results may be summarized as follows: 

There exist two critical chemical potentials fic(e) and fic(e), ^c\e) < fie \e). 

For jj, < ijLc\e), the solution of the variational problem ()2.36|) is oo^^. For the state 

a + = a_ = a;, = 0, (4.1) 

i.e. the two a — ± Bose gases behave like two mean field Bose gases with no BEC and they 
do not interact with the external 6-boson laser field, in which there is no condensation either. 

/q\ /r\\ ————— 

For n > fie '(e), the solution of the variational problem is ' . By virtue of (j3.13|) for this 
state we have 

< a + < a_ and a b ^ 0, (4.2) 

i.e. there is BEC for the two a — ± Bose gases and for the external boson laser field 
(superradiance). Moreover, for e > the condensation of the a — ± bosons is enhanced by 
the presence of this laser field: one gets it even for dimensions v — 1, 2, because po(— s) < oo 
for v > 1. We interpret this quantum state as that of equilibrium BEC superradiance 

When p^\e) < /4 (e), for jJi (e) < fi < ^c\e), the solution of the variational problem is 
ujp ' . For this state we have 

a_ 7^ and a + = a b = 0, (4.3) 



i.e. there is BEC only for the a = — Bose gas. 

(a) The simplest case to consider is when e > x . In this case ^c\e) = fii(e) and ^c\s) = 
/^(e). Thus the states are uj^\ ujp and uo^ as [i increases. This means that as we increase 
/j, we observe three stages: the system goes from no BEC, to BEC for only the a = — bosons, 
and then to BEC for both a — ± boson species and for the laser field (superradiance). 

(b) When < e < x we have to consider two subcases: fi±(e) < yU < /-^(s) and /i < fii(e). 

In the first subcase H\(e) = ^c\e) < /jtf\s) < At 2 (e). Otherwise the situation is as in (a). 

For /i < fii{e), < t^c\e) < /42(e) but we did not determine if ^c\e) < fi¥\s). Thus 

we do not know if the intermediate phase with BEC for only a = — bosons is present. 

(c) If e = 0, then ^c\e) = l^i\s) and the intermediate phase with BEC for only a = — 
bosons is not present. 
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